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We investigate giant resonances of spherical nuclei on the basis of the Argonne V18 potential after uni- 
tary transformation within the Similarity Renormalization Group or the Unitary Correlation Operator Method 
supplemented by a phenomenological three-body contact interaction. Such Hamiltonians can provide a good 
description of ground-state energies and radii within Hartree-Fock plus low-order many-body perturbation the- 
ory. The standard Random Phase Approximation is applied here to calculate the isoscalar monopole, isovector 
dipole, and isoscalar quadrupole excitation modes of the 40 Ca, 90 Zr, and 208 Pb nuclei. Thanks to the inclu- 
sion of the three-nucleon interaction and despite the minimal optimization effort, a reasonable agreement with 
experimental centroid energies of all three modes has been achieved. 

PACS numbers: 21.30.Fe,21.45.Ff,21.60.Jz 
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I. INTRODUCTION 



The most consistent starting point for nuclear structure the- 
ory are nuclear Hamiltonians derived from quantum chromo- 
dynamics (QCD) in the framework for chiral effective field 
theory containing two- and three-nucleon interactions 01 Ho- 
using these interactions we can employ unitary transforma- 
tions, e.g. the Similarity Renormalization Group (SRG) or 
the Unitary Correlation Operator Method (UCOM), to pre- 
diagonalize the Hamiltonian and to improve the convergence 
behavior of various many-body approaches. Recently, this ap- 
proach was applied successfully to light and medium-mass 
nuclei in the context of the No-Core Shell Model 0, 0] and 
in Coupled-Cluster Theory and related methods ]5h7|] . 

The computational effort, however, limits the applicability 
of general three-nucleon interactions in the unitary transfor- 
mation as well as in the application in many-body methods. 
Furthermore, to provide an appropriate starting point for the 
investigation of collective excitations in the framework of the 
Random Phase Approximation (RPA) the interaction has to 
reproduce experimental ground-state radii reasonably well — 
this has not yet been achieved on the basis of chiral two- plus 
three-nucleon interactions beyond the lightest nuclei. As a 
preparatory step towards the full inclusion of chiral two- and 
three-nucleon interactions, we follow a more pragmatic ap- 
proach by using the unitarily transformed Argonne V18 po- 
tential [8] supplemented by a phenomenological three-body 
contact interaction. This allows us to investigate ground- 
state properties as well as collective excitations throughout the 
whole nuclear mass range up to 208 Pb. In a previous paper the 
influence of phenomenological three-nucleon interactions on 
the description of ground-state nuclear properties was inves- 
tigated [21 and a good simultaneous description of ground- 
state energies and radii was achieved. This was not possible 



with the pure two-body UCOM interaction employed in ear- 
lier studies ifioll . We now examine whether dynamical prop- 
erties, such as collective modes, show a similar quantitative 
improvement compared to previous work [fTTTl . Therefore, in 
this Brief Report we apply the standard RPA to study collec- 
tive excitations of closed-shell nuclei. We will show that the 
results obtained with the unitarily transformed Argonne VI 8 
supplemented by a phenomenological three-body contact in- 
teraction agree within 20% with the results for traditional phe- 
nomenological potentials like the Gogny D1S interaction IU2I1 . 
The characterization of these hybrid Hamiltonians in standard 
applications is mandatory to provide a well-defined footing 
for predictive calculations like the study of low-energy dipole 
transitions [13], and investigations in the framework of quasi- 
particle RPA (QRPA) fll . 

After a presentation of the formalism in Sec. [II] we dis- 
cuss briefly some ground-state properties of closed-shell nu- 
clei across the whole nuclear chart on the basis of Hartree- 
Fock and many-body perturbation theory in Sec. [Ill] In Sec. 
|IV| we study the isoscalar monopole (ISM), isovector dipole 
(IVD), and isoscalar quadrupole (ISQ) excitation modes of 
three chosen nuclei, 40 Ca, 90 Zr, and 208 Pb. 



II. FORMALISM 

The unitary correlation operator method (UCOM) and the 
similarity renormalization group (SRG) provide two differ- 
ent approaches for the generation of soft phase-shift equiv- 
alent two-body interactions. These two methods have already 
been discussed extensively (see B15II and refs. therein). Based 
on the Argonne V18 potential we will apply four different 
classes of unitarily transformed interactions in the follow- 
ing: UCOM(SRG), S-UCOM(SRG), SRG, and S-SRG which 
were introduced in Jgl]. These transformed two-body interac- 
tions are supplemented by a simple three-body contact inter- 
action 
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with variable strength C^. The m-scheme matrix elements of 
the contact interaction can be evaluated on-the-fly, which is 
of great computational advantage. As single-particle basis the 
eigenstates of the harmonic oscillator are employed. For cal- 
culations beyond the mean-field level we have to introduce a 
regularization of the contact interaction which is achieved via 
a restriction of the total oscillator energy of the three-particle 



state: {2n\ + h) + (2« 2 + h) + (2«3 + h) < £3 



where n 



and I are the principal and the angular momentum quantum 
numbers of the harmonic oscillator states, respectively |]9|]. 



III. GROUND-STATE PROPERTIES 

For a first characterization of the four different two- plus 
three-body interactions, the Hartree-Fock (HF) approxima- 
tion and many-body perturbation theory are used to calculate 
ground-state energies and charge radii for selected closed - 
shell nuclei from 4 He to 208 Pb. The formal inclusion of the 
three-body contact interaction in these two methods has been 
discussed in detail in Ref. J3]. The single-particle basis 
is truncated with respect to the principal oscillator quantum 
number e — 2n + 1, which is restricted to e < e max = 14 with an 
additional constraint for the orbital angular momentum quan- 
tum number I < / max = 10. The oscillator parameter ano is 
chosen for each nucleus separately such that the experimental 
charge radius is reproduced by a shell-model Slater determi- 
nant built from harmonic oscillator single-particle states. We 
have observed that for these values of «ho and for e max = 14 
the HF ground-state energies deviate by less than \%c from 
their minimum with respect to «ho and that the mean square 
radius is stable against variations of «ho- Thus the HF ground- 
state energy and mean-square radius are well converged with 
respect to e max . 

Figure Q] shows the ground-state energies per nucleon and 
charge radii for closed-shell nuclei across the whole nuclear 
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FIG. 1: Ground-state energies per nucleon and charge radii of se- 
lected closed-shell nuclei resulting from HF calculations (filled sym- 
bols) and MBPT (open symbols) based on the UCOM(SRG) (•), S- 
UCOM(SRG) (■), SRG (♦), and S-SRG (a) interactions using the 
optimal parameter sets (cf. Tab.|IJ. The bars indicate the experimen- 
tal values (13, [H]. 





a [fm 4 ] 


C 3N [GeVfm 6 ] 


UCOM(SRG) 


0.16 


1.6 


S-UCOM(SRG) 


0.16 


2.2 


SRG 


0.10 


4.3 


S-SRG 


0.10 


2.0 



TABLE I: Optimal parameter sets for the different two- plus three- 
body interactions. 



chart. The three-body cut-off parameter is set to Ei max = 20 
for all four interactions (cf. Q). The strength of the three- 
body contact interaction is used to adjust the charge radii to 
the experimental values. As observed in Fig. [T]the strength 
C3N can be chosen such that the radii are well reproduced 
across the whole mass range by all four interactions. This is 
a remarkable result considering the simplistic structure of the 
three-body interaction. The optimal values for the three-body 
contact strengths are summarized in Table U together with the 
corresponding flow parameters. 

The HF ground-state energies reproduce the systematics of 
the experimental values except for an almost constant shift, 
which is due to the missing effects of long-range correlations 
that cannot be described at the mean-field level. The influ- 
ence of long-range correlations can be taken into account via 
second-order many -body perturbation theory lHoj . Here, we 
apply perturbation theory only to the two-body part of the in- 
teraction (cf. H]). The effect on charge radii is negligible, 
but for the ground-state energies the inclusion of the second- 
order perturbative corrections leads to a substantial improve- 
ment (open symbols in Fig.Q]). The agreement with the exper- 
imental data is not yet perfect, i.e. the differences vary from 
0.2 to 2.4 MeV per nucleon, but one has to keep in mind that 
the energy corrections are not yet fully converged with respect 
to the single-particle model space size. Furthermore, we only 
consider the second-order estimate and have no information 
about the influence of the third and higher orders lH8l U9ll . 

In case of the UCOM(SRG) interaction one observes that 
the energy corrections for some nuclei, mainly the heavier 
nickel isotopes and 100 Sn, do not follow the general trend. 
The origin of this behavior is found in the corresponding HF 
single-particle spectra, where some level spacings are col- 
lapsed leading to divergent contributions to the perturbative 
corrections. Therefore, these data points are not shown here. 



IV. GIANT RESONANCES 

For the investigation of collective excitations we apply the 
standard Random Phase Approximation (RPA) (cf. e.g. lUlll ). 
In the framework of the HF approximation and the standard 
RPA the three-body contact interaction is equivalent to the 
density-dependent two-body interaction ||2fl-|2jl 



Vnn[£>] 



C3N 
6 



(2) 



which is used for the implementation of the RPA for compu- 
tational reasons. The only difference to the HF plus MBPT 
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calculations discussed in the previous section is the absence 
of the three-body cut-off £3 msx in the RPA. The HF calcula- 
tions providing the basis for RPA are also performed using the 
density-dependent form of the three-body interaction without 
the cut-off. As the HF energies are independent of Ej, max the 
chosen implementation has no influence on the RPA results. 
Although we use density-dependent two-body interaction in 
RPA we will still refer to it as three-body contact interaction. 

In the following, we investigate three excitation modes, 
namely isoscalar monopole (ISM), isovector dipole (IVD), 
and isoscalar quadrupole (ISQ) excitations, of 40 Ca, 90 Zr, and 
208 Pb. The different two- plus three-body interactions intro- 
duced in the previous sections are considered using the pa- 
rameter sets listed in TableQ] The single-particle basis is again 
truncated at <? max = 14. 

As a first benchmark we consider the exhaustion of the 
classical sum rules. We validate our implementation by us- 
ing the Gogny D1S interaction, whose momentum dependent 
terms are of zero range: the ISM and ISQ classical sum rules 
are then fulfilled within 1 % or better. In case of the trans- 
formed AVI 8 interactions, the exhaustion of the classical ISM 
sum rule lies between 90.5% for 40 Ca calculated with the 
UCOM(SRG) interaction and 98.4% for 208 Pb calculated with 
the S-SRG interaction. The exhaustion of the ISQ sum rule 
lies between 98.4% for 208 Pb calculated with the SRG inter- 
action and 102.9% for 40 Ca calculated with the UCOM(SRG) 
interaction. The deviations from the classical IS sum rules are 
mostly a consequence of the non-localities of the finite-range 
interactions used ll23ll . When including the three -body inter- 
action we use a larger flow parameter in order to compensate 
for the additional repulsion. The larger flow parameter gener- 
ates stronger non-localities in the transformed two-body inter- 
action, which in turn influence the exhaustion of the classical 
sum rules. It is expected that the Thomas-Reiche-Kuhn sum 
rule for the IVD mode is significantly enhanced due to the 
non-localities of the applied interactions. Indeed, the smallest 
percentage that we found for the Thomas-Reiche-Kuhn sum 
rule was 168.3%, for 40 Ca based on the SRG interaction. Fi- 
nally, we have found that the energy of the spurious dipole 
state never exceeds 20keV in the present cases. 

In Figure [2] we compare the HF and RPA response func- 
tions for all three excitation modes calculated with the S- 
UCOM(SRG) interaction for 208 Pb. The response functions 
are obtained via a convolution of the calculated discrete 
strength distribution with a Lorentzian function with a width 
of 2 MeV. The upper row shows the response functions ob- 
tained with the pure two-body S-UCOM(SRG) interaction us- 
ing the flow parameter a = 0.04 fm 4 while the response func- 
tions depicted in the lower row were obtained with the two- 
body S-UCOM(SRG) interaction plus the three-body con- 
tact interaction with a = 0.16 fm 4 and = 2.2GeVfm 6 
(cf. Tab. H]i. The HF response is spread wide in case of both 
isoscalar modes while it is rather compressed for the isovector 
excitation. In comparison, the RPA response is compressed 
significantly and shifted to lower excitation energies in the 
isoscalar channels which leads to strongly collective excita- 
tion modes, the giant resonances. In case of the ISQ mode 
one observes the excitation of a low-lying 2 + state in addition 





(a) 


(b) 


(c) 


(d) 


Exp. 


(e) 


ISM 


40 Ca 


19.81 


18.57 


18.91 


17.92 


19.18 


21.06 




90 Zr 


16.73 


16.40 


17.37 


15.81 


17.81 


17.53 




208p b 


12.88 


12.93 


14.41 


12.35 


14.18 


13.19 


IVD 


40 Ca 


25.86 


23.22 


20.34 


22.61 


21.9 


22.70 




90 Zr 


22.20 


20.79 


18.46 


20.27 


17.9 


19.82 




208p b 


17.88 


16.96 


15.72 


16.44 


13.6 


16.21 


ISQ 


40 Ca 


22.54 


20.06 


19.02 


19.92 


17.8 


17.69 




90 Zr 


18.49 


17.45 


16.75 


17.38 


14.2 


12.70 




208p b 


14.07 


13.66 


13.43 


13.54 


10.9 


9.36 



TABLE II: Centroid energies in MeV obtained for the (a) 
UCOM(SRG), (b) S-UCOM(SRG), (c) SRG, and (d) S-SRG inter- 
actions, compared to experimental values j24l - l29ll and the (e) Gogny 
D1S interaction. 



to the giant resonance. The RPA response of the IVD excita- 
tion is shifted to higher excitation energies with respect to the 
HF response, i.e. the residual interaction is attractive in the IS 
channels and repulsive in the IV channel. 

Comparing the HF response functions obtained with the 
pure two-body interaction with those resulting from the two- 
plus three-body interaction reveals that the inclusion of the 
three-body interaction leads to a compression of the response 
for all considered excitation modes. This compression can 
be understood by considering the HF single-particle spectra 
which are spread wide when calculated with a pure two-body 
interaction. The repulsion of the three-body interaction in- 
creases the level density and thus leads to a compression of 
the HF response. 

The RPA response calculated with the two- plus three-body 
interaction is concentrated in a narrower resonance structure 
compared to the response functions obtained with the pure 
two-body interaction for all three excitation modes. Further- 
more, the centroids are shifted to lower energies: the ISM 
centroid is shifted by 1 MeV while the IVD and ISQ centroids 
are moved by 3 MeV, respectively. 

Figure [3] summarizes the response functions that were ob- 
tained with the four different two- plus three-body interac- 
tions for 40 Ca, 90 Zr, and 208 Pb including centroid energies 
extracted from experiment. All four interactions yield com- 
parable results with only minor differences and are in reason- 
able agreement with the experimental centroids. The ener- 
gies of the giant dipole and quadrupole resonances tend to be 
overestimated, but much less so than with the pure two-body 
UCOM interaction ITTll . The systematic nature of these devi- 
ations hint at the role of higher-order configurations beyond 
RPA (H. 

The results obtained with the S-UCOM(SRG) and the S- 
SRG interactions are very similar in all aspects confirming the 
similarities between these two interactions that were already 
observed earlier 

The centroid energies of the calculated response functions 
are compared to the experimental values and to calculations 
based on the Gogny D1S interaction in Table HH In most of 
the 36 cases the centroids are within 20% of the experimental 
value. The SRG interaction performs particularly well while 
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FIG. 2: Comparison of HF ( ) and RPA ( ) response functions calculated with the pure two-body S-UCOM(SRG) interaction (upper 

row) and the two-body S-UCOM(SRG) plus three-body contact interaction (lower row) for 208 Pb. The ISM, IVD, and ISQ response functions 
are given in units of fm 4 / MeV, e 2 fm 2 / MeV, and e 1 fm 4 / MeV, respectively. 




the largest deviations are observed for the UCOM(SRG). Our 
results are similar in quality with existing results using phe- 
nomenological energy density functionals, e.g. the Gogny 
D1S interaction listed in Table UJ We note that the good over- 
all agreement of the response with experiment is obtained de- 
spite the systematic deviation of the HF ground-state energy 
(without MBPT corrections) from the experimental binding 
energy. The same remarkable observation was made in [fl4ll . 
It is well known that an excellent description of nuclear ener- 
gies by phenomenological energy density functionals based 
on HF does not guarantee a good description of collective 
states within RPA. For example, they usually require a high 
nucleon effective mass, which leads to a strong underestima- 
tion of the IVD energy. Our present and previous results lUill 



consistently confirm that the HF or HF-Bogolyubov ground- 
state energy has little to do with the dynamical behavior of 
the system under study, as described by RPA or quasi-particle 
RPA. Instead it is more important to obtain a good description 
of the radii and the ground-state energies beyond the mean- 
field level, e.g. in many-body perturbation theory. 



V. CONCLUSIONS 

We have investigated the ISM, IVD, and ISQ giant reso- 
nances using unitarily transformed two-nucleon interactions 
supplemented by a simple phenomenological three-body con- 
tact interaction. For the three considered nuclei, 40 Ca, 90 Zr, 



5 



and - 08 Pb we achieved a reasonable agreement with experi- 
mental centroid energies, where the SRG interaction yielded 
the smallest deviations. We have shown that a reasonable de- 
scription of charge radii is important to obtain a good starting 
point for RPA calculations. 

Following this benchmark, the hybrid interactions em- 
ployed here can be applied for the quantitative description of 
other collective phenomena in closed and open-shell nuclei. 
Several follow-up studies are in progress. At the same time, 
our results represent an intermediate step towards a consis- 
tent inclusion of chiral two- plus three-nucleon interactions 
for the study of collective excitations. Apart from computa- 
tional challenges, the systematic reproduction of ground-state 



radii will be a key issue when working with these Hamiltoni- 
ans. 
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